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Abstract 

A statistical theory is developed for the stochastic Burgers equation in the inviscid limit. Master equations 
for the probability density functions of velocity, velocity difference and velocity gradient are derived. No 
closure assumptions are made. Instead closure is achieved through a dimension reduction process, namely 
the unclosed terms are expressed in terms of statistical quantities for the singular structures of the velocity 
field, here the shocks. Master equations for the environment of the shocks are further expressed in terms of 
the statistics of singular structures on the shocks, namely the points of shock generation and collisions. The 
scaling laws of the structure functions are derived through the analysis of the master equations. Rigorous 
bounds on the decay of the tail probabilities for the velocity gradient are obtained using realizability 
constraints. We also establish that the probability density function Q{^) of the velocity gradient decays 
as 1^1"^''^ as ^ — CX3. 
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1 Introduction 



Consider the randomly forced Burgers equation: 



(1.1) 



Ut + UUx = VUxx + /, 



where f{x,t) is a zero-mean, Gaussian, statistically homogeneous, and white-in-time ran- 
dom process with covariance 



and B{x) is a smooth function. In the language of stochastic differential equation (1.1) 
must be interpreted as 



We will be interested in the statistical behavior of the stationary states (invariant mea- 
sures) of (1.1) if they exist, or the transient states with possibly random initial data. 
There are two main reasons for considering this problem. The first is that (1.1) and its 
multi-dimensional version are among the simplest nonlinear models in non-equilibrium sta- 
tistical mechanics. As such they serve as qualitative models for a wide variety of problems 
including charge density waves [19], vortex lines in high temperature super-conductors [5], 
dislocations in disordered solids and kinetic roughening of interfaces in epitaxial growth 
[29], formation of large-scale structures in the universe [36, 40], etc. The connection be- 
tween these problems and (1.1) can be understood as follows. Consider an elastic string in 
a random potential V{x, s). The string is assumed to be directed in the sense that there is 
a time-like direction, assumed to be s, such that the configuration of the string is a graph 
over the s-axis. Let Z{x,t) be the partition function for the configurations of the string 
in the interval < s <t, pinned at position x at time t: 



where {-{wit) = x) denotes the expectation over all Brownian paths w{-) such that w{t) = 
X, P = 1/kT, k is the Boltzmann constant, T is the temperature. The free energy, 
'p{x,t) = In Z{x,t) then satisfies 



{f{x,t)f{y,s)) = 2B{x-y)6it-s), 



du = (—uux + i^Uxx)dt + dW{x, t), 



where W{x,t) is a Wiener process with covariance 



{W{x, t)W{y, s)) = 2B{x - y) min(t, s). 




where u = kT. This is the well-known Kardar-Parisi-Zhang equation [27]. In one dimen- 
sion if we let u = ipx^ f = ^) we obtain (1.1). 
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The second reason for studying (1.1) is in some sense a technical one. For a long 
time, (1.1) has served as the benchmark for field-theoretic techniques such as the direct 
interaction approximation or the renormahzation group methods, developed for solving 
the problem of hydrodynamic turbulence. This role of (1.1) is made more evident by the 
recent flourish of activities introducing fairly sophisticated techniques in field theory to 
hydrodynamics [33, 8, 25]. In this context (1.1) is often referred to as Burgers turbulence. 
Since the phenomenology of the so-called Burgers turbulence is far simpler than that of 
real turbulence, one hopes that exact results can be obtained which can then be used 
to benchmark the methods. However so far our experience has proved otherwise: The 
problem of Burgers turbulence is complicated enough that a wide variety of predictions 
have been made as a consequence of the wide variety of techniques used [1, 2, 3, 4, 6, 7, 
9, 10, 11, 12, 13, 15, 22, 23, 24, 25, 26, 28, 33, 35, 37, 38, 40]. 

The main purpose of the present paper is to clarify this situation and to obtain exact 
results that are expected for Burgers turbulence. Along the way we will also develop some 
technical aspects that we believe will be useful for other problems. The main issues that 
interest us are the scaling of the structure functions and the asymptotic behavior of the 
probability density functions (PDF) in the inviscid limit. The former is well-understood 
heuristically but we will derive the results from self-consistent asymptotics on the master 
equation. The latter is at the moment very controversial, and we hope to settle the 
controversy by deriving exact results on the asymptotic behavior for the PDFs. 

From a technical point of view we insist on working with the master equation and 
making no closure assumptions. The unclosed term is expressed in terms of the statistics 
of singular dissipative structures of the field, here the shocks. We then derive a master 
equation for the statistics of the environment of the shocks by relating them to the singular 
structures on the shocks, namely the points of shock creation and collisions. These are 
then amenable to local analysis. In this way we achieve closure through dimension reduc- 
tion. We then extract information on the asymptotic behavior of PDFs using realizability 
constraints and self-consistent asymptotics. We certainly hope that this philosophy will 
be useful for other problems. 

One main issue that will be addressed in this paper is the behavior of the PDF of the 
velocity gradient. Assuming statistical homogeneity, let (5'^(^,t) be the PDF of ^ = Ux- 

satisfies 

Qt = iQ"" + i^^Ql^ + BiQ'i^ - i^ii^xxlOQn^, 

where Bi = —Bxx{0), {^xx\0 is the average of ^xx conditional on ^. This equation is 
unclosed since the explicit form of the last term, representing the effect of the dissipation, 
is unknown. We are interested in Q'^ at the inviscid limit: 

Q{^,t) = limQ''{^,t). 
— I/— >o 

In order to derive an equation for Q, one needs to evaluate 



F{^,t) = -limu{{^xx\OQnt 
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This is where the difficulty arises. 

Remembering that Ux{x,t) = lim./\x^o{u{x + Ax,t) — u{x,t))/ Ax, the procedure out- 
Hned above pertains to the process of first taking the limit as Ax 0, then the limit as 
— > 0. It is natural to consider also the other situation when the limit ^ is taken 
first. In this case the limiting form of (1.1), (notice that we now write the nonlinear term 
in a conservative form), 

(1-2) nt + Uu\ = f, 

has to be interpreted in a weak sense by requiring 

(1.3) jj dxdt ^uipt + \u^'^x + /y'} = 0, 

for all compactly supported smooth functions ip. The solutions u satisfying (1.3) are called 

weak solutions. In order to ensure a well-defined dynamics for (1.2), i.e. existence and 
uniqueness of solutions with given initial data, an additional entropy condition has to be 
imposed on weak solutions. This amounts to requiring 

(1.4) u{x+,t) < u{x-,t), 

for all {x,t). The entropy condition (1.4) is the effect of the viscous term in the inviscid 
limit. There is a huge mathematical literature on (1.2) for the deterministic case. Standard 
references arc [30, 31, 39]. The random case was studied recently in the paper by E, 
Khanin, Mazel and Sinai [14]. 

The first step in the present paper is to derive master equations for single and multi- 
point statistics of u satisfying (1.2). In particular, we derive an equation for the PDF of 
ri{x,y,t) = {u{x + y) — u{x,t))/y, Q^{ri,x,t). We are interested in 

Qi^,t) = ]imQ\C,x,t). 
x— >o 

One natural question is whether 

(1.5) Q = Q. 

We will present very strong argument that (1.5) holds for generic initial data and for the 
type of forces described after (1.1). Some of our results also apply to Q for the more 
general case when it is possibly different from Q. 

The issue now reduces to the evaluation or approximation of F{S^,t). Several different 
proposals have been made, each leads at statistical steady state {Qt = 0) to an asymptotic 
expression of the form 



Q(0 



C-ier" as C^-oo, 

C+^/?e-«V(3Bi) as C^+oo, 
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but with a variety of values for the exponents a and (3 (here the C±'s are constants). By 

invoking the operator product expansion, Polyakov [33] suggested that F = aQ + bS,Q, 

with a = and b = —1/2. This leads to a = 5/2 and /3 = 1/2. Boldyrev [6, 7] considered 

the same closure with — 1 < 6 < 0, which gives 2 < a < 3 and /? = 1 + 5. Based on heuristic 

arguments, Bouchaud and Mezard [9] introduced a Langevin equation for the local slope 

of the velocity, which gives 2 < a < 3, /? = 0. The instanton analysis [4, 18, 25] predicts 

the right tail of Q without giving a precise value for /?, and it does not give any specific 

prediction for the left tail. E et al. [13] made a geometrical evaluation of the effect of F, 

based on the observation that large negative gradients are generated near shock creation. 

Their analysis gives a rigorous upper-bound for a: a < 7/2. In [13], it was claimed that 

this bound is actually reached, i.e., a = 7/2. Finally Gotoh and Kraichnan [24] argued 

1/3 

that the viscous term is negligible to leading order for large |^|, i.e. F « for |^| 3> , 
giving rise to a = 3 and /3 = 1. 

In this paper, we will first give a more detailed proof of the bound 



by studying the master equation for the environment of the shocks. 

An important reason behind the success of this program is that local behavior near 
the most singular structures (here the points of shock creation) is understood. Points of 
shock creation were isolated in [14] as a mechanism for obtaining large negative values of 
Ux and resulted in the prediction that a = 7/2. Prom this point of view, the present paper 
provides the missing step establishing the fact that points of shock creation provide the 
leading order contribution to the left tail of Q. 

For the most part, our working assumption will be that solutions of (1.2) are piecewise 
smooth. Shock path are smooth except at the points of collision. In particular, shocks are 
created at zero amplitude and the shock strength adds up at collision. For transient states 
these statements follow from standard results for the deterministic case under appropriate 
conditions for the initial data [32]. For stationary states, these results are of the type 
established in [14] under a non-degeneracy condition for the forcing. 

Before ending this introduction, we make some remarks aboTit notations and nomen- 
clature. In analogy with fluid mechanics u will be referred to as the velocity field. We will 
denote the multi-point PDFs of u{-, t) as Z'^{ui,xi, . . . , Un,Xn, t), i.e. 



a > 3, 



announced in [15]. We will then show that 



a = 7/2, 



Prob(ai < u{xi,t) < 6i, . . . , a„ < u{xn, t) < 6„) 




dui--- / dUn Z''{ui,Xi,...,Un,Xn,t). 



The superscript refers to the viscous case, > 0. In the inviscid limit we will denote the 
multi-point PDFs of u{-,t) by Z{ui,xi . . . ,Un,Xn,t). Statistical stationary values will be 
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denoted by the subscript oo, e.g. or Z^. We reserve the special notations R''{u, x, t) = 
Z'^{u,x,t) and R{u,x,t) = Z{u,x,t) for the one-point PDF of u{x,t). Q'^{C,t) denotes 
the PDF of (^{x,t), and its inviscid Umit is Q{^,t). Statistical symmetries, or equality in 
law, will be denoted by '= . Two of such symmetries will be repeatedly used. The first is 
statistical homogeneity, 

(1-6) u{x, t) '= u{x + y, t) for all y. 

(1.6) holds if for all measurable, real-valued functional ^{u{-,t)) we have 

= mu{- + y,t))) for ally. 

Note that from (1.1) and the assumptions on /, it follows that (1.6) holds for all times if 
uo{x) '= uo{x + y), where fxo(-) = u{-,0) is the initial velocity field. Also, if a statistical 
steady state exists and is unique, then it satisfies (1.6). The second symmetry will be 
referred to as statistical parity invariance, and is related to the invariance of (1.1) under 
the transformation x — >■ —x, u — > —u. This implies that 

(1-7) u{x,t)^= —u{—x,t), 

or 

Mn(;t))) = {^i-ui-i-),t))). 

(1.7) holds for all times if uo{x) '= —uo{—x). (1.7) is also satisfied at statistical steady 
state. 

Even though some of our results are formulated in terms of lemmas and theorems, the 
emphasize here is on the ability to calculate things rather than rigor. We will assume that 
the initial distribution Zq{ui,xi, . . . ,Un, Xn) = Z{ui,xi, . . . , n„, 0) is concentrated on 
D{I), the space of functions defined on / C M admitting only jumps discontinuities: 

D{I) = {'u(x, t), ^(xit, f) exists for allx, with u{x+,t) < u{x-,t)}. 

Note that the existence of stationary states is only established when / is a finite interval 
on M [14]. We will also assume that the initial velocity field uo{-) = u{-,0) is independent 
of ly and statistically independent of /(•)• 

Finally, unspecified integration ranges are meant to be (— oo,-|-oo). 

2 Velocity and velocity diflFerences 

Let u'^{x,t, f,uo) be the solution of (1.1) with forcing /, initial data uq. It follows from 
standard results (for the deterministic case, see [30, 31]) that for fixed t 

u''{-,t,f,uo)^u{-,t,f,uo) asu^Q, 
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P X Po-almost surely, in L^{I), where I is the interval on which (1.1) is studied. Hence 



weakly. For the statistically stationary states, it was established in [14] that, weakly, we 
have 

Therefore to study the inviscid limit of statistical quantities of u, such as the PDFs of u 
and Su{x, y, t) = u{x + y,t) — u{y, t), it is enough to consider (1.2). 

It is useful, however, to write this equation in a modified form more convenient for 
calculations. To this end, let 

u±{x, t) = u{x±, t). 
For any function g{u) = g{u{x,t)), we define two average functions 

bWlA = U9{u+) + g{u-)) 

(2.1) 



b(«)]B = [ dp 9{u- + f3{u+ - u-)). 
Jo 



It is shown in [41] that (1.2) is equivalent to 

(2.2) ut + [u]j^Ux = /, 
or, in the language of stochastic difi^erential equations, 

(2.3) du = -[u]j^U3:dt + dW{x,t). 

(2.1) assigns an unambiguous meaning to the quantity [nj^n^; when the solutions of (2.2) 
develop shocks. Moreover the following chain and product rules hold 

(2.4) g3;{u) = [gu{u)]j^Ux, 



(2.5) {g{u)h{u))^ = [g{u)]^h^{u) + [h{u)]^g^{u). 

Similar rules apply for derivatives in t. As an example of these rules, we have that the 
integral of the term [ti]^^! across a shock located at a; = y is given by (using [u]^ = [u]^) 
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u{x, t) = t) + u-{x, t)), 

s{x,t) = — U-{x,t). 

In terms of (u, s), the two averages in (2.1) are 

rl/2 

[g{u)\^ = / df3g{u + lis). 

7-1/2 

We will denote by {yj} the set of shock positions at time t. Hence the set 
quantifies the shocks at time t. 

2.1 Master equation for the inviscid Burgers equation 

We now turn to R{u,x,t), the one-point PDF of the solution of (2.2). We have 

Theorem 2.1 R satisfies 

(2.6) Rt = -uR^ - j du' K{u - u')Rx{u', x, t) + BqRuu + G, 

where Bq = B(0), K{u) = {H{u)—H{—u))/2, H{-) is the Heaviside function and G{u, x,t) 
is given by 

/ rO ru-s/2 \ 

(2.7) G{u,x,t) = \ ds du {u — u)Q{u,s,x,t) ] . 



— oo Ju+s/2 

Here g{u, s, x, t) is defined such that g{u, s, x, t)dudsdx gives the average number of shocks 
in + dx) with u{y,t) G [u,u + du) and s{y,t) G + ds), where y E [x,x + dx) is 
the shock location. 

Remark 1. G can be referred to as a dissipative anomaly. It can be written more 
explicitly as 

1 /-o 

G{u, x,t) = - ds s{q{u — s/2, s, x, t) + q{u + s/2, s, x, t)) 

(2-8) ^1/2 .0 

— / dp ds sg{u + Ps, s,x,t). 

J-l/2 J-oo 



STATISTICAL THEORY FOR BURGERS EQUATION 10 
Remark 2. g{u,s,x,t) can be equivalently defined as 

where 



(2.6), (2.7) do not require statistical homogeneity. For homogeneous situations, these 
equations simphfy. Since Rx = (2.6) reduces to 

(2.9) Rt = BqRuu + G, 

where G{u,x,t) = G{u,t). Since the shock characteristics are independent of its location, 
we have 

g{u,s,x,t) = pS{u,s,t) 

where p = p{t) is the number density of shocks, S{u,s,t) is the PDF of {u{yo,t),s{yo,t)), 
conditional on the property that yo is a shock position (because of the statistical homo- 
geneity, yo is a dummy variable). Thus, (2.7) reduces to 

(rO ru-s/2 \ 

ds du {u-u)S{u,s,t) . 

J-oo Ju+s/2 J ^ 

At statistical stationary state, i?t = in (2.9), and one readily verifies from this 
equation that 



(2.11) Roo(u) = —— / ds du \ — -{u-uf\ Sooiu, s), 

IHq J-oo Ju+s/2 \ 4 J 

Clearly, i^oo > since S(X){u, > 0. In addition, by direct computation we obtain 

J du i?oo(«) = 

Thus, from the requirement that R^o be normalized to unity, we get 

(2.12) Bo = -^{s^). 
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Proof of Theorem 2.1: Let 9{\,x,t) = e"*'*'"^^'*) . {9) is the characteristic function 
of u{x, t) and R is given by 

R{u,x,t)= I ^e'^-{9{X,x,t)). 

Using Ito calculus, 

dW{x, t)dWiy, t) = 2B{x - y)dt, 

it follows from (2.3), (2.4), (2.5) that 

d9 = {-i\[u\^u^[9]^ - X^Bo9)dt - i\9dW{x,t). 

Thus 

{9)t = iX{[u],u,[9l)-X'Bo{9). 

Note that the terms involving the force contain 9, not This is because 9 = except 

when there is a shock, and this is a set of zero probability. Of course this argument does 
not apply for the convection term since Ux is infinite at shocks. To average the convective 
term «A['u]^nx[0]B we use 

Ox = -i>^Ux[9]j, 

to get 

iX{[u]^Ux[9l) = -{[u]Jx) = -{u9)^ + {ux[9]^). 

For convenience we write this equation as 

iX{[u],uMn) = -i{9U + i^-\0)x + {uM, - [9]^)). 

Combining these expressions gives: 

(2-13) {9)t = -i{9)xx + i\-\e)x - X'Bo{9) + G, 

where G{X,x,t) is given by 

G{X,x,t) = {uMA-[OW)- 

To proceed with the evaluation of G, note that the only contributions to this term are 
from the shocks, since [9]j^ = [0]^ except at shocks. Let {yj} denote the positions of the 
shocks at time t. Since Ux{yj, t) = s{yj,t)S{x — yj) at the shocks, G can be understood as 

^(A,^,*) 



"1/2 

d 

-1/2 
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where sj = s{yj,t). Using (2.1), 
G{X,x,t) 

= (Y1 - yj>~'^^' U^'i/^ + e-^^^^/2 - 2 r 

where Uj = u{yj,t). Using g{u,s,x,t) this average is 

= J duj^ ds _^^-iXs/2 _2 j^'\(3 ^-'^1^'^ q{u, s,x,t), 

Going back to the variable u, we get (2.8), hence (2.6). 

I 

Remark. Under the assumption of ergodicity with respect to spatial translations, an 
alternative derivation of (2.10) is to go back to (2.14) and use the equivalence between 
ensemble average and spatial average. Then 

G(A,t) 

TV 1 ^ 
i=i 

where N is the number of shocks in the interval [—L,L\. Using the ergodicity again it 
follows that the sum is equal to 

= p j du y° ds|e-^^^(^e^^^/2 + e-^^^/2 _ 2 e"^^^"^ S{u,s,t), 

where we used limL_>oo -/V/(2L) = p. Going back to the variable u, we get (2.10). 



2.2 An alternative derivation 

The derivation of (2.6) (or (2.9) for statistically homogeneous situations) given above is 
rigorous, but rather unintuitive. In fact, the effect of the viscous term is buried in the 
definition of the two averages in (2.1), and it is not clear at all at this stage whether (2.6) 
arises in the limit of the equation for as — > 0. Assuming statistical homogeneity, 
recall that R^ satisfies (see (A. 4) in the Appendix) 

Rt=B^Rl^-u{{u,,\u)R^)^. 
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Here we give another, less rigorous but more intuitive, derivation of (2.9) and (2.10) by 
working with the equation for , and calculating directly the limit of the viscous term 
as z/ — > 0. In particular, we will compute explicitly that the dissipative anomaly (2.10) is 

given by 

(2.15) G{u,t) = -\hnu{{uxx\u)R)^. 

This will give strong support to the claim that R = lim,y^o R'^ satisfies (2.9). 

Assuming spatial ergodicity, the average of the dissipative term can be expressed as 

v{uxx\u)R'' = v{uxx{x,t)5{u - u{x,t))) 

(2.16) 1 rL 



1 

= v lim — / dx Uxx{x,t)5{u — uix^t)). 

L—*oo 2L J— I, 



Clearly, in the limit as v ^ only small intervals around the shocks will contribute to 
the integral. In these intervals, boundary layer analysis can be used to obtain an accurate 
approximation of u{x,t). 

The basic idea is to split u into the sum of an inner solution near the shock and 
an outer solution away from the shock, and using systematic matched asymptotics to 
construct uniform approximation of u (for details see, e.g., [21]). For the outer solution, 
we look for an approximation in the form of a series in u: 

u = u°"* = uo + uui + 0{y^). 

Then uq satisfies 

UQt + UqUQx = f, 

i.e. Burgers equation without the dissipation term. In order to deal with the inner solution 
around the shock, let y = y(t) be the position of a shock, and define the stretched variable 
z = {x — y)/^ and let 

u'''{x,t) = V (^^^ + S,?j , 

where 6 is a perturbation of the shock position to be determined later. Then, v satisfies 
(2.17) uvt + {v-u + v^)vz = Vzz + J^f, 

where u = dy/dt, 7 = dS/dt and, to 0{u'^), vf can be evaluated at x = y (and can thus 
be considered function of t only). 

We study (2.17) by regular perturbation analysis. We look for a solution in the form 



V = Vq + UVl + 0{l'^). 
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The boundary condition for this equation arises from the matching condition with = 
uo + vui + 0(z/^): 

hm Vo = hm uq = u± - , 

z— >±oo x^y 2 

where s = s{t) is the shock strength. It is understood that for smah ly matching takes 
place for small values of |a; — y| and large values of \z\ = \x — yl/u. This gives 

s / sz 
vq = u — - tanh f — 

These results show that, to 0{i'), (2.16) can be estimated as 



J 



where Qj is a layer centered at yj with width S> 0{u). Going to the stretched variable 
z = {x — y)/v, and taking the hmit as L — >■ oo, we get 

//•O /'+00 
du / ds S{u,s;t) / dz vozziz:'t)^{''J' ~ vo{z,t)), 
J —oo J —oo 

where S{u,s;t) is the PDF of {u{yo,t), s{yQ,t)) conditional on yo being a shock location. 
The z-integral can be evaluated exactly using 

dzvozz = dvo^^^ = dvo{vo - u), 

VOz 

where we used the equation [vq — u)voz = vq^z- This leads to 

//•O /•u—s/2 
du / ds S{u,s;t) / dvo (vo — u)d{u — vq). 
J —oo Ju+s/2 



Hence, 

ds 

-oo Ju+s/2 

Using this expression in (2.15), we get (2.7). 



/O /•u—s/2 
ds / du {u — u)S{u, s;t). 
-oo Ju+s/2 
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2.3 Computing the anomalies 

Here we derive equations for the moments of R. For simplicity, consider a statistically 
homogeneous case, and assume that u{x,t) '= —u{—x,t). Then S{u,s,t) = S{—u,s,t), 
G{u,t) = G{—u,t), R{u,t) = R{—u,t). This means that all moments of odd orders of 
R are zero. For the moments of even orders, we get from (2.9) the following equations 
(n G No): 

^(u^") = 2n(2n - l)So(tx2-2) + h2n, 



where /i2n is the anomaly term 
h2n = J du u^^'G 



^ ' "m - sf'^{u + ns)) - ((2« + sf^iu - ns 



22"(2n + 1) 

(^2n+i = by parity). An alternative definition of h^n is 



h2n = 2n lim iy{u^''~\^^) = -2n(2n - 1) lim iy{u^"~'^ul). 



This gives for instance 



/i2 = -2 1im^n2) = ^(/). 



At statistical steady state, it gives 

(^,2n^ = %^(((2tx - s)2"+2(tx + (n + l)s)) - {{2u + s)2"+2(u - (n + l)s 



Bo 

where Cn = n\/2'^^'^(n + 3)!. These expressions can also be obtained from (2.11). In 
particular, for n = 2, we obtain again (2.12). 

2.4 Multi-point PDF 

We now turn to Z{ui,xi, . . . ,Un,Xn,t), the multi-point PDF of u. We have 
Theorem 2.2 Z satisfies 

n n 
Zt = — ^ UpZxp + ^ B{Xp — Xq)Zu^Uq 

n „ 

(2.18) ~S / du' K{Up-u')Z^^{ui,Xi,...,u',Xp,...,Un,Xn,t) 

p=l'' 

n 

+ ^ ^ G(Up, Xp, 112, X2 ■ ■ ■ , Up— I, Xp—i, Up-^-i, Xp-^.\, ■ ■ ■ , Un, Xji, t), 
p=l 
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where K{u) = {H{u) — H{—u))/2, H(-) is the Heaviside function and G is given by 



gives the average number of shocks in [xi,xi + dx) with u{y,t) G [u,u + du), s{y,t) G 
[s,s + ds) where y G [xi,xi + dx) is a shock location, and with u{x2,t) G [u2,U2 + du2),- ■ ■ , 

u{Xn, t) G [Un, Un + dUn) ■ 

We will omit the proof of Theorem 2.2 since it is a straightforward generalization of 
the one given for Theorem 2.1. 

For the two-point PDF, Z{ui,xi,U2,X2,t), (2.18) becomes 



G{ui,Xi, . . .,Un,Xn,t) 



(2.19) 




Here g{u, s,xi,U2,X2 ■ ■ ■ , Un, Xn, t) is defined such that 



q{u, s, xi,U2, X2, ■ ■ ■ , Ufi, Xnt)dudsdu2 ■ ■ ■ dundx, 



(2.20) 




+ 2B{xi - X2)Zu^u2 
+G{ui,Xi,U2,X2,t) + G{U2,X2,Ui,Xi,t). 



Assuming statistical homogeneity, (2.20) can be simplified. First, 



Z{ui,y, U2,x + y,t) = Z{ui,0, U2,x, t) = Z{ui,U2, x, t). 



Secondly, 



q{u, s, y,u,x + y, t) = pT{u, s, u, x, t) 



where p is the number density of shocks, and T{u, s,u,x,t) is the PDF of 



{u{yo,t),s{yo,t),u{yo + x,t)), 



STATISTICAL THEORY FOR BURGERS EQUATION 17 

conditional on yo being a shock position. Thus for statistically homogeneous situations, 
(2.20) reduces to the following equation for Z{ui,U2,x,t) 

Zt = -{U2 - Ui)Zx 

— J du' K(u2 — u')Zx{ui,u' ,x,t) 
(^•^^) + j du' K{ui-u')Zx{u',U2,x,t) 

+BoZuiUl + BoZu2U2 + 2-B(x)Zu^u2 

+G{U1,U2,X, t) + G{U2,UI, -X, t), 

where 

G{ui,U2,X,t) 

(2.22) / /"O j-u-i-s/2 _ _ _ \ 

= P\ ds s du {ui — u)T{u, s, U2, x,t) \ 

\J —oo Jui+s/2 J u\ 

2.5 Velocity difference and structure functions 

Assume statistical homogeneity, and let Z^{Su,x,t) be the PDF of the velocity difference 
Su{x, y, t) = u{x + y,t) — u{y, t). Z^{w, x, t) is related to Z{ui,U2, x, t) by 



Z^{w, x,t) = J du Z ^ — — , 



w 

u + —,x,t 



It then follows immediately from (2.21) and (2.22) that 
Corollary 2.3 Z^ satisfies 

^2 = -^4 -2jdw' H{w' - w)Zi{w\ x, t) 

+2{Bo-B{x))Zi^ + G\w,x,t), 

where K{w) = {H{w) — H{—w))/2, H{-) is the Heaviside function and 



(2.24) 



G^{w, x,t) = J du G ^u — ^,u+ —,x,t 

+ J du G (u + ^,u-^,-x,t^ 



Remark. G can be put into a form which is more convenient for the calculations. Let 

Su+{x,yo,t) = u{yo + \x\,t) - u+{yo,t), 
6u-{x,yo,t) =u-{yo,t) - u{yo - \x\,t), 
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and let U±{s,Su±,x,t) be the PDFs of {s{yo,t),Su±{x,yo,t)) conditional on y being a 
shock position. Then, can be expressed as 

G\w, X, t) = G%{w, X, t) + Gt{w, X, t), 

where 

. p fO 

G^{w,x,t) = - / ds s{U±{s,sgn{x)w — s,x,t) + U±{s,sgD.{x)w,x,t)) 

^ J —oo 

— pi ds s dp U±{s,sgn{x)w — Ps,x,t). 

J-oo Jo 

We now consider some consequences of (2.23). We have the following two theorems: 
Theorem 2.4 In the limit as x ^0, 

Z\w,x,t) = d{w) - \x\{p6{w) + iOS^w) - pS{w,t)) +o{x), 

(2.25) ^ 1 /w \ 

Z%w, X, t) = {l- p\x\)-Q \x\pSiw, t) + o{x), 

where 5^{w) = d5(w)/dw, o{x) must be interpreted in the sense of weak convergence. Here 
Q{^,t) is the PDF of^{x,t), the regular part of the velocity gradient, i.e. 

Ux{x,t) = ^{x,t) +Y,siyj,t)S{y - yj), 

3 

and S{s,t) is the PDF of s{yo,t) conditional on yo being a shock location. 
Theorem 2.5 (Structure function scahng) Let 

(2.26) (|(5n|") = J dw \w\''Z\w,x,t). 
In the limit as x — > 0, 

r |xr(|er) + o(a;«) if 0<a<l, 

(2.27) {\6un = \x\m) + p{\s\)) + oix) ^/ a = 1, 

[ Ixlpdsl") + o{x) if 1< a. 

In terms of the moments, (2.27) is (n G No) 

.2 28) {6u'n = \x\p{s'n + o{x), 
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The first moment satisfies {6u) = for all {x,t). This is a consequence of statistical 
homogeneity and it is readily verified since multiplying (2.23) by w and integrating gives 
in {6u)t = 0. 

An equation for Q{^,t) will be derived in Sec. 3. To interpret (2.25) note that can 
be decomposed into 

Z^(w,x,t) = pnsix,t)Z^ (w, x,t\no shock) 

+ {1 - Pnsix,t))Z^{w,x,t\shock), 

where Pns{x,t) is the probability that there is no shock in [y,y + x), Z^{w,x,t\no shock) 
is the PDF of 6u{x,y,t) conditional on the property that there is no shock in [y,y + a;), 
Z^{w,x, 1 1 shock) is the PDF oi6u{x,y, t) conditional on the property that there is at least 
one shock in [y, y + x). Since, by definition of p we have 

Pns = 1 - P\x\ +o{x), 

(2.25) states that 

Z^{w, X, t\no shock) = (1 — p\x\)—Q ^—,t^ + o{x), 
Z\w, x,t\shock) = S{w,t)+o{l). 

This is consistent with the picture that Su{x, y, t) = x^{y, t) + o{x) if there is no shock in 
[y, y + x), and w{x, y, t) = s{yo, t) + o(l) if yo G [y,x + y) is a shock position. 

Proof of Theorem 2.4: Wc will consider the case x > 0. The case x < can be 
treated similarly. Note first that, in the limit as x ^ 0, we have 

xZ\xi,x,t)^Q{i,t), 

weakly. We postpone the proof of this fact until Sec. 3. It implies that 

Z\w,x,t) = 5{w) + o{l), 

weakly. Define 

Aiw^t) = \im. x~^ ( Z^ (w,x,t) — 5(w)) = Mm. Z^(w,x,t). 

Taking the limit as x — > in the equation for Z^ , it follows that A satisfies 

(2.29) = -wA -2 j dw' H{w - w')A{w', x, t) + B{w, t), 

where we used lima,^o(-So — B{x)) =0 and we defined 

B{w,t) = lim G^{w,x,t). 
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To evaluate B note that as x ^ 

6u±{x,yo,t) 0, 

almost surely. This implies that, as a; — 0, 

U±{s,w,x,t) — > S{s,t)5(w), 

where S{s,t) is the PDF of s{yo,t) conditional on yo being a shock location. Hence, from 
the expression for G^, 

/w 
dw'S{w't) - 2pH(w), 
-oo 

where H{-) is the Heaviside function and we used S{s,t) = for s > since s{yo,t) < 0. 
Inserting this expression in (2.29), the solution of this equation is 

A{w,t) = -5{w) + p{s)6'^{w) + pS{w,t). 

Here 5^{'w) = d5{w)/dw and we used the identity w6^{w) = —6{w). Using (3.4), p{s) = 
— (^), this can be restated as 

A{w,t) = -S{w) - {OS\w)+pS{w,t). 

Hence, combining the above results, we have 

Z\w,x,t) = S{w)-x{S{w) + {^)S^{w) - pS{w,t)) + o{x), 

weakly. This establishes the first equation in (2.25) for a; > 0. Reorganizing this expression 
as 

(2.30) Z^{w,x,t) = {1 - px){S{w) - x{OS^{w)) + xpS{w,t) + o{x), 

weakly, and using the identity 

5{w) - x{^)5\w) = +oix), 

establishes the second equation in (2.25) for x > 0. 

I 

Proof of Theorem 2.5: Wc will prove (2.27) directly for moments of integer order 
higher than one. For other values of a (2.27) follows from (2.25) and the fact that the tails 
are controlled by higher order moments. Note first that for a > 

lim(|,5nr) = 0, 

X—^0 
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since 6u{x, y,t) ^ almost surely as x ^ 0. Now, multiply (2.23) by (n G N, n > 
integrate and take the limit as x — Oib. The result is 



where 



Note that 



a^= lim X ^((5n")= lim 

x—*0± >0± 

6^ = lim / dw w'^G\w,x,t). 
I dw w"G\w,x,t) = ^(s((5u+ +sgn(x)s)") + ^{s6ul) 

J Zi Zi 

- p f^dp {s{Su+ + Psga{x)s)'') 

Jo 

+ ^{s{6u.+sgnix)sy') + ^{s5u'l) 

- p Cdp (s((5u_ +/3sgn(x)sf ). 

JO 

Since 5u±{x, yo,t) — > almost surely as x — > 0, we have 

bt = (±1)>(.'^+^) (l - ^) . 
Inserting this expression in the equation for gives 

a±+i = ^hm x-i(5u"+i) = (±l)>(."+i). 

Thus 

{Su^'') = |x|p(s2") + o(x), = + o{x). 

This proves (2.28). 

We now prove (2.27) for < a < 1. The proof for other values of a is similar. Let 

fiw, X, t) = {l- p\x\)^Q (^^, tj + \x\pS{w, t), 
g\w,x,t) = S{w) - \x\{pS{w) + {05\w) - pS{w,t)), 
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and write for M > 

= / dw {wl'^iZ^ -f) + I dw \w\''{Z^ - f). 

J\w\<M J\w\>M 

The first term at the right hand-side is o{x) because of (2.25). To estimate the second 
term, note that for M large enough 

dw \w\''Z^ < I dw w'^Z^ 

\w\>M J\w\>M 



( dw w^{Z^ - + / dw w^g' 

J J\w\>M 

= o{x) + \ApI dww^Siw,t) 

J\w\>M 

= o{x) + 0{x) / dw w'^S{w, t). 

J\w\>M 



I dw \w\-f = |xr (1 - \x\p) I di ierg(e, o 

J\w\>M J\€,\>M/x 

+ \x\p / dw \w\°-S{w,t) 

Jw>M 

= o(x") + 0{x) I dw Iwl^Siw, t). 

Jw>M 

Since M can be made arbitrarily large, we get 

J dw Iwl^iZ^ - f) < o(x") + 6mO{x), 

where (5jv^ — >■ as M — >■ +00. Noting that 

kr(i^r) if < a < 1, 



J dw Iwl'^f 
we obtain (2.27) for < a < 1 



\x\m) + p{\s\)) if a = l. 



3 Velocity gradient 

We now turn to the study of the PDF of the velocity gradient for the solution of (2.2). 
Since the solutions typically contain discontinuities, it is already an issue whether the 
PDF for the velocity gradient is well-defined. Heuristically, this is so since u only fails to 
be diffcrcntiable at no more than countably many points. We will therefore be concerned 
with the PDF of the regular part of the gradient. 
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3.1 Master equation for the PDF of the velocity gradient 

We focus on the statistically homogeneous case and derive an equation for Q{^,t), the 
PDF of ^{x,t), defined as the regular part of the velocity gradient, i.e. 

Ux{x,t) = ^{x,t) + YjS{yj)S{x - yj), 

3 

where i{-,t) G L^{I)- We will prove that 
Theorem 3.1 Q satisfies 

(3.1) Qt = iQ + {eQ)^ + B^Q^^ + F{i,t), 
where Bi = —B^xi^) o,nd 

(3.2) Fi^,t) = p r dssV{s,C,t). 

J —oo 

Here V{s,^,t) = {V+{s,^,t) + V.{s,^,t))/2, V±{s,^±,t) are the PDFs of 

{s{yo,t),^±{yo,t) = Ux{yo±,t)), 

conditional on the property that yo is a shock position. 

The consequences of (3.1) will be studied in Sec. 3.3. Note that if u{x,t) '= —u{—x,t), 
then 

iC+{x,t),s{x,t)) = {^^{-x,t),s{-x,t)), 

andV+{s,^,t) = V_{s,(,t) = Vis,C,t). 

Proof of Theorem 3.1: Let Q^{ri,x,t) be the PDF of r]{x,y,t) = {u{x + y,t) - 
u{y,t))/x. is related to Z and by 

du Z iu — + — . 

From (2.23) it follows that satisfies 

Ql = VQ' + {v'Q\ + Bi{x)Qf,^ - XT^Qi 
^^■^^ -2x J drl H{r] - r]')Qi{r]', x, t) + F\r), x, t), 

where Bi{x) = 2{Bo — B{x))/x'^ and 

F\r], X, t) = Ffiri, x, t) + F^{ri, x, t) + ^^(r?, x, t), 
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with 



Ff{ri,x,t) = p f ds sV^{s,ri,x,t), 
F2 {rj, x,t) = p J ds sV^ ^s, 7] — -j — j-, X, t 



\x\ 



Fi{ri,x,t) = -2p y° ds s d(3 (^s,ri - j^^,x,t 

Here V\s,rf,x,t) = {Vf{s,rf,x,t) + Vl{s,r],x,t))/2, Vl{s,r),x,t) are the PDFs of 

{s{yo,t),r]±{yo,x,t)), 

with ri^{yQ,x,t) = ±{u{yo it \x\,t) — u±{yQ,t))/\x\, conditional on the property that yo is 
a shock position. 
Define 

Q{^,t) = limQ\^,x,t). 
x— »o 

It is easy to see that 

J dri F^{ri,x,t)=p{s), 

drj F^{r],x,t) = —2p drj ds sV^ {s,r],x,t) = —2p{s), 



■00 



consistent with the fact that 



and, hence, 

d_ 
dt 



j dv F\r],x,t)=0, 
j dj] Q^T], X, t) = 0, j dj] TjQ^r], x, t) = 0. 



In the hmit as x ^ 0, 77+ and r/_ converge respectively to the gradient of the velocity at 
the left and right sides of the shock, and, pointwise in 77, we have 

lim Fi{n,x,t) = lim Fi{n,x,t) = 0, 
lim F^{ri,x,t)=Firi,t), 
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with F given by (3.2). Therefore, a standard argument with test functions apphed to (3.3) 
shows that in the hmit as x ^ 0, converges weakly to Q, solution of (3.1). 

I 

Remark. F2 and arc examples of what we will call "ghost terms", i.e. terms that 
have finite total moments, but in the limit converge pointwise to zero. Due to the ghost 
terms, is not clear at the moment that Q satisfies / Q{^,t) = 1. This will be established 
as a consequence of Lemma 3.2. 

3.2 Alternative limiting processes 

Using BV-calculus, one works at = and access the statistics of the velocity gradient 
by taking x — ^ in {u{x + y,t) — u{y, t))/x. This procedure gives an equation for 

Q(e,t) = limQ^(e,x,t), 

a;— >0 

where Q^(^, x, t) is the PDF of (u(x + y, t) — u{y, t))/x. In this section, we revert the order 
of the limits: we take x ^ first, working at finite ly, and then let ly ^ 0. As in Sec. 2.2, 
this is done using boundary layer analysis and matched asymptotics. In this way, we will 
obtain an equation for 

Q(^,t) = limg^(e,t), 
— I/— >o 

which will turn out to be identical to the one for Q. To the extend that boundary layer 
analysis can be justified, this strongly suggests that the limits x ^ 0, — > commute. 

For statistically homogeneous situations, recall that Q'^{^,t) satisfies (see (A. 5) in the 
Appendix) 

Qt = ^Q" + i^^Qli + BiQl^ - 

The average of the dissipative term can be expressed as 

1 

As in Sec. 2.2, we will evaluate this integral using the approximation for ^xx provided by 
the boundary layer analysis. However, this analysis has to be developed further in order 
to evaluate the dissipative term for the velocity gradient. Recall that = u^^ = v^s^jv + 
v\z + 0{y\ Inside the shock, the 0(1) contribution of vi^ is clearly negligible compared 
to the 0(i^~^) contribution of vq^/v. However, the contribution of is important at the 
border of the shock because vq^ decays exponentially fast there. 

To evaluate vi^ at the shock boundaries z ±00, consider the equation for vi that we 
get from (2.17): 
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VI = Cie./'o'^^'(^o(^')-") 

+ f^dz' {C2 + {vo{z') - nh + dz"vot{z") - el>" (M^'-), 

where Ci, C2 are constants. They, as well as 7, have to be determined by matching with 

^out _ _|_ _|_ 0{v^) [21] . We will not dwell on this problem since Ci, C2, 7 do not enter 
the expression for vi^ as z —i- ±00. Indeed, using the expression for vq, direct computation 
shows that vi reduces asymptotically to 



1 / du ds „ , 



and 



1 f du ds „\ 
2 / ds du s^7 



s2 



Thus 



2du Ids 2/ , 

hrn vi^ = T--j7- --jT± — = ^±, 

0-»±oo s dt s dt s 

where the last equality is simply a definition of Note that these can be reorganized to 
give 

ds s ,^ ^ . du s ,^ X ^ 

In the limit as ^ these are the equations of motion along the shock. 
Using these results, to 0{v), (2.16) can be estimated as 

^{ixx\^)Q'' = j dsdud^+d^- S{u,s,C+,^-,t) 

X / dx e^:c{x,m^ - r{x,t)), 

Jn 

where Q is a layer centered at y with width » 5(?2, s, i) is the PDF of 

{u{yo,t), s{yo,t),^+{yo,t),^-{yo,t)), conditional on j/o being a shock location (^_|_ and ^_ 
have to be included for reasons to be made clear later). 

We now go to the stretched coordinate z = {x — y)/^, and use the result of Sec. 2.2, 
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namely, 

^- = Z.-l^;0_,+t;i^ + O(l/). 

To 0(z/), both terms must be included. However the contribution of vi^ is important only 
at the border of the shock where vq^ falls exponentially fast, and can be neglected inside 
the shock. Thus, « ^ = i'~^vqz + and we have 

/+0O 
-oo 



To perform the z integral, we change the integration variable to ^' = ^ — ^_ for z < and 
to ^' = - C+ for z>0. For z < 

and we have 
Similarly, for z > 

{ivo-u)i^-^+))z = ^zz, 

K -u) = - (Is' + 2Ke - ^+)) = I + 0{i^), 

and we have 



This leads to 



= p [ dsdud^+d^^ S{u, s, ^+,^-,t) f d^ ^ - ^' - 

+ p I dsdud^+d^. s{u,s,^+,^.,t) f di' '-5{i - e' - 
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Letting u —>■ and integrating gives 

ds s 

-oo 

where we used the consistency constraint 



/O roo 
dss / d^' Vis,^',t), 
-oo J£ 



Thus 



j dud^^S{u,s,$,+ ,$,_,t) = V±{s,^±,t). 



3.3 Consequences of the master equation 

First we observe 
Lemma 3.2 

(3.4) {^)+p{s)=0, or Jd^m^,t) = -p{s). 

Proof: Denote by u'^ the solution of (1.1) for finite u. Then (u^) = 0. Let (p{x) be 
a compactly supported smooth function. We have 



= J dx ip{ul) = l^j dx ipul^ = -l^j dx ifxvyy 
In the limit as i/ ^ 0, ^ u, the solution of (2.2). Thus 

= —(^J dx <PxU^ ■ 
Denote by yj the location of the shocks. We can write 

f dx ipxU = - XI / dx ipux - ip{yj)s{yj,t). 



'Vj 

Averaging this result, we get 



dx ^{{e)+p{s)) = Q, 

for all compactly supported smooth functions (p. Hence (3.4). 
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Notice that for finite u or x, f d^^Q^{^,t) = f d^^Q^{^,x,t) = 0. In the language 
of Kraichnan [28], (3.4) represents a flow of probability of ^ from the smooth part of 
the velocity field to the shocks. It reflects the fact that no matter how small v is, the 
dissipation range has a finite effect on inertial range statistics. 

As a consequence of (3.1) and (3.4), we have 



i.e. the normalization of Q is preserved. For the initial data we are interested in, this 
implies 



In particular J ^F(^,t) = p((s^+) + (s^_))/2 exists and is finite. 

Proof: Since p{s) and its derivative are obviously finite, the existence of the integral 
/ d£, ^F{^, t) = p{{s^+) + (s^_))/2 follows from the argument below and a standard cut-off 
argument on large values of 





Lemma 3.3 



(3.5) 



Using 



u±t + ^^±^± = /, 



where ^±{x,t) = Ux{x±,t), and dyj/dt = u{yj,t), we have 



—u+{yj,t) = -^C+{yj,t)+u+t{yj,t) 

= u{yj,t)^+{yj,t) - u+{yj,t)^+{yj, t) + /{yj, t) 
= -2^{yj^'t)^+{yj^'t) + f{yj,t). 



Similarly 



—u-{yj,t) = -s{yj,t)C-{yj,t) + f{yj,t). 



These equations can be reorganized into 



^u{yj,t) = ~s{yj,t){^+{yj,t) - i-{yj,t)) + f{yj,t). 
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and 

Consider now 

(3.6) p{s) = (y2^i^^*)^(^ -yju = -%■)/• 

j j 

Using the equation for s{yj,t) and dyj/dt = u{yj,t), we have 

|(p(^)) = -f((^e+) + (^U) 



-( s{yj,t)uiyj,t)5{x - yj) 
j 

+contribution from shock creation and colhsion. 



The second term at the right hand-side is zero by homogeneity. Thus to obtain (3.5) it 
remains to prove that the contribution of shock creation and collision vanish. To under- 
stand how this term arises and why it vanishes, assume a shock is created at position yi 
at time ti. Then the sum under the average in (3.6) involves a term like 

Ti = s{yi,t)6{x-yi)H{t-h). 
where -ff (•) is the Heaviside function. Time differentiation gives 

= ^s{yi,t)5{x - yi)H{t - h) - s{yi, t)u{yi,t)5x{x - yi) 
+ s{yi,t)5{x - yi)5{t - ti). 

The last term accounts for shock creation. Since the shock amplitude is zero at creation, 

s{yi,t)6{x - yi)6{t - ti) = s{yi,ti)6{x - yi)6{t - ti) = 0. 

This means that shock creation makes no contribution to the time derivative of (3.6). 
Consider now the collision events. Assume at time ti the shocks located at y2 and 2/3 
merge into one shock located at yi. Obviously yi(ti) = 1/2(^1) = Vsiti)- Such an event 
contributes to the sum under the average in (3.6) by a term like 



T2 = s{yi,t)6{x-yi)H{t-h) 

+ {s{y2, t)5(x - ys) + s(y3, t)5(x - y^))H{ti - t). 
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Time differentiation gives 

'^ = ^^s{y,m^-y,)H{t-t,) 

+ ^s(2/2, t)5{x - y2)H{ti -t) + ^s{y3, t)6ix - yz)H{ti - t) 

- s{yi,t)u{yi,t)5.^{x ~ yi)H{t - h) 

- {s{y2,t)u{y2,t)6^{x - y2) + s{y^,t)u{y^,t)5^{x - ys))H{ti - t) 
+ is{yi,t)5{x - yi) - s{y2, t)6{x - ^2) - s(2/3, t)6{x - y3))6{t - h). 

The term involving S(t — ti) arises from shock colUsion. Since yi(ii) = 2/2(^1) = 2/3(^1) 
shock amphtudes add up at collision, 

^lim s{yi{h +t),h+ t)S{x - yiih + t)) 

= ^'^}^_^i^(.y2{'ti - - - y2{ti - 1)) 

+s{y3{ti -t),ti- t)5{x - y^ih - t))), 

the term in the equation for T2 involving S{t — ti) vanishes. This means shock collision 
makes no contribution to the time derivative of (3.6). Hence (3.5). I 



Corollary 3.4 We have 
(3.7) 



hm |C|'Q(e,0 = 0, 



i.e. Q decays faster than |^| ^ as ^ — > — 00. 

Proof: Taking the first moment of (3.1) leads to 



d_ 
dt' 

Using (3.4) this equation can be written as 

d 



+00 



+ I ((s^+) + (s^_)) 



dt 



{p{s)) 



P 



eQ -'-{{si+) + {si.)). 



Using (3.5) we obtain that the boundary term must be zero. Since S^^Q has different sign 
for large negative and positive values of ^, one must have separately lim^_^_oo \i\^Q = 
and lim^^+00 \^?Q = 0. Hence (3.7). I 

Remark. Lemma 3.3 uses the fact that shocks are created at zero amplitude and shock 
strengths add up at collision, which holds in the case of large-scale smooth forcing [14]. It 
is possible to construct pathological situations, such as the unforced case with piecewise 
linear initial data [34], in which case shocks are created at finite amplitude. In this case 
Lemma 3.3 is changed to 
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Lemma 3.5 

(3-8) |(p(.)) = -^((s^+) + (.U)-^c, 

where 

Dc = -Cri(si) - (72 (S2) > 0. 

Here ai and cj2 o,re respectively the space-time number density of shock creation and col- 
lision points, (si) < is the average shock amplitude at creation, and {32) < is the 
average gain of amplitude at shock collision. In this case 

(3.9) g(e,t)~Z),|e|-3 as ^^-00. 

Lemma 3.5 follows from a direct adaptation of the proof of Lemma 3.3. 

3.4 Realizability and asymptotics for the statistical steady state 

We now turn to the study of (3.1) at steady state 

(3.10) = ^Q+ieQ)^ + BiQ^^ + F{0, 

where F{^) = lim^^oo F{^, t)- We shall prove that 

Theorem 3.6 The realizability constraint Q G L^(M), Q >0 implies 

(3.11) lim rVF(e) = 0, 
^ ' 4— ++00 

where A = —^^/2>Bi. Assuming (3.11), the only positive solution of (3.10) can he expressed 



as 



(3.12) Q^iO = ^l^ < i'FiO - ^ < ^""'Gia 



where 

Furthermore 

(3.13) Qoo(0 
where 



i>l J-OO 




^ d^^'Fi^') as e ^ -00, 

as ^ — > +00, 
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Under the conditions of Theorem 3.6 we also have 

Lemma 3.7 Qoo can be expressed as 

? / J^' ^/-2„A' / 7/-// i:llT?lcll\ 



(3.14) Qoo(e) = -v/ ^c'rv/ de"me 

for ^ < 0, and 

(3.15) QUi) = C+^e-^ - ^ /^V r'e^' /^"^i^" ^"^^(^'0, 
for e > 0. 

Remark 1. ^F{^) is integrable as a consequence of Lemma 3.3. In particular (3.5) at 
steady state gives 



Remark 2. C+ is finite if (3.11) holds. To see this, note that because of the factor e^, 
problem may arise at ^ = +oo only. (3.11) implies that we can write 

i^(0 = e-V(0, 

with /(^) = o(^2) as C ^ +00. Using Jd^^F = 0, we write G{^) as 

^ f + OO 



c r+oo 

t r+oo 

= e-V(e)-^/ <e'e-^7(0. 

-Dl 

Using e~^' = —Bi^'~^{e~^')^' , after integration by parts we have 

G(0 = / d^' (r V(C'))«'e-^ . 



Since /(O = o{e) as ^ ^ +oo, {^'fiOh = and 

/•+00 , / r+co ,\ / \ 

GiO = d^ o{l)e-^' = a e-^'j = a (r'e"^) , 

where at the last step we used e~^ = 0((^~^e~^)^). Thus 

e^G(0 = o(r') as ^ ^ +oo, 
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which imphes that C+ is finite. Similarly, we can show that the last integral in (3.15) is 
finite iff (3.11) holds. Indeed, using F{^) = o{^^e~^) as ^ ^ +oo we have 



= o(e'e-^'), 

where we used ^^e~^ = 0((^e~^)^). Thus 

(3.16) r'e^ / df fFif) = o(r ') as ^ ^ +oo, 

which is the necessary and sufficient condition for the last integral in (3.15) to be finite. 

Proof of Lemma 3.7: To show (3.14), we start from the explicit expression for the 
integral involving G in (3.12) 

-Dl J-OO 

-Ol J —oo J^i J —oo J —oo 

and integrate by parts the second integral using e^' = Bi^'~^{e~^')^i . The result is 



-/ <e^'G(^') 

>1 J-oo 



Bi 

Inserting this expression in (3.12) gives (3.14). 

To show (3.15), we use / d£,^F = and write (3.12) as 

1 r+°° fp~^ /■+°° ./ 

Qoo(6 = C+^e-^ - — / d^' ^FiO + ^ / < G(C')- 
-Dl r>i 

Using J d^^F = we write explicitly the integral involving G in this expression as 



Bi 



e^'F(eo - ^ di'e'^'i'j^^ dee Fin 
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and integrate by parts the second integral using e^' = Bi^'~^(e~^')^f: 



Inserting this expression in the last expression for Qoo gives (3.15) 



Bi Ji Bi Jf' 



Proof of Theorem 3.6: The general solution of (3.10) is Q = Q0D + C1Q1 + C2Q2, 
where Ci and C2 are constants, Qi and Q2 are two linearly independent solutions of the 
homogeneous equation associated with (3.10). Two such solutions are 

(3.17) Qi(0=ee-^, 

(3.18) Q,(^) = i_^ d^'^'e^\ 

-Dl J-QO 

For ^ < 0, using e^' = 5i^'~^(e^')^/, after integration by parts Q2 can be written as 

(3.19) Q2iO = -^e-^ f d^' r'e^' 

J —00 



We now show that realizability requires Ci = C2 = 0. First, one readily checks that 
lim^_,._oo Qoo = lim^->-oo Q2 = 0, while Qi grows unbounded as ^ — > —00. Hence in order 
that Q be integrable we must set Ci = and the general solution of (3.10) is 

Q{0=C2Q2iO + Qoo{0- 

We evaluate this solution asymptotically as ^ ^ ibcxD. Consider Q2 for large negative ^ 
first. Using e^ = Bi^'~^(e^ )^/, after integration by parts we write (3.19) as 

J —00 

Since ^~^e^ = 0{{^~'^e^)^) as ^ — —00, the integral in this expression is of the order 

ee-^ d^' r'e^' = ee-^ d^'0{{r'e\) 
J —00 J— 00 

= oUe-^ f d^'iC'e\) 
^ J —00 ' 
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Thus 

Q2{o = Slier' + o(r') as e ^ -oo, 

Consider now Qoo for large negative ^. Using e^' = Bi^'~'^(e^')^i , after integration by 
parts we write (3.14) as 

J— oo 

J—oo ^ J —oo 



As e — — OO 

J—oo J—oo ? J—oo 

^ — oo 

where we used the estimate above. Thus 

Qoo(6 = ICr' i'F{^') + o(r') as C ^ -oo. 

J—oo 

Combining these expressions, we have 



Q(e) = C2i?i|a-3 + ia-3/^ di'^'FiO 

J—oo 



(3.20) 

+ 0(C2r') + o(r') as C ^ -oo. 



Consider now Q2 for large positive ^. Write (3.18) as 

Q2(o = 1 - ^ e'e^' - ^ /V e'e^', 

-Dl J -00 -Dl Jf* 



where > is arbitrary but fixed. Using e^' = Bi^' ^(e^')^/, after integration by parts 
of the second integral we get 
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Integrating by parts again gives 

Q2(6 = -BiC' - 4Siee-^ /V r'e^' + ©(^e-^). 

Since ^~^e^ = 0((^~^e^)^) as ^ ^ +cxd, the remaining integral can be estimated as above, 
yielding 

Q2(o = -5ir'+o(r')- 

Considering Q^o for large positive ^, we must distinguish two cases. If (3.11) does not 
hold, then in (3.12) we decompose 

^ [' < e^'G(e') = ^ e^'GiO + ^ /V e^'G(e') 

= ^ /V e^'G(0 + 0(^e-^). 

and write 



e ^1 



Qooio = -B- 1 < i'm') - 1 e^'G(e') + o(e 

The second integral in this expression is explicitly 

-Dl 



-Dl J -CO 



Bi 

The integration by parts of the second integral using e^' = Bi^'~'^{e~^')^i gives 



- / di' e^'Gii') 
1 ^e* 



1 /•€ 



4/ '^^'^'^(o 

-Dl J-oo 



- /^c^e' e'-'e^' de eF{o + 0(ee-^) 

1 J^i, J-oo 



B 

Thus 



Qoo{0 = 1%' C'-'e^' di" fFia + 0(^e-^^ 

-Dl J^^, J-oo 
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Integrating by parts again using e^' = Bi^' ^(e^')^/ gives 



J —CO 



The last integral can be estimate as for the large negative ^. Using / d£,^F = 0, this leads 
to 



Ar+oo 

and, hence, 



r+oo 



r+cc 



(3.21) 

+0(C2r')+o(r') as ^^+oo. 

In contrast, if (3.11) holds, then we can use (3.15). From (3.16) it follows that the integral 
term in this expression is o(^e~^) as ^ — ^ +oo. Thus, 

Q(e) = C+ee-^ + o(ee-^), 

and 

(3 22) = -C2Bir' + C+^e-^ 

+0(C2r^) + 0(^6-"^) as e ^ +00. 

In (3.20), (3.21), (3.22), if the C2 term at the right-hand side is non-zero, then it will 
dominate the second term. However, since the C2 term has opposite sign when ^ — > ±00, 
it must be zero, i.e., we must set C2 = 0. This proves that Q = Qoo- Furthermore, since 
the F term at the right-hand side of (3.21) is negative (recall that F < 0), this solution 
must be rejected in order that Q be non-negative. Thus (3.11) must hold, and we get 
(3.13). 



Remark. If the assumptions of Lemma 3.3 do not apply, i.e. if shocks are not created 
at zero amplitude or shock strength does not add up at collision, then from (3.8) at steady 
state we have 



/ d^ mi) = -De. 
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This implies that (3.21) has to be changed to 
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r+oo 

+o(C2r') + o(r') as e^+oo. 



This sohition cannot be dismissed as being unreaHzablc since the second term at the right 
hand-side can balance the first one. In particular, for C2 = D(./B\ we obtain from (3.20) 
Q(^) = Dc\^\-'^ + o(^-3) as ^ -cx), consistent with (3.9). 

3.5 Statistics for the environment of the shocks 

We now turn to the statistics for the environment of the shocks. For simplicity we will 
focus on statistically homogeneous situations such that u{x,t) '= —u{—x,t). Define 



conditional on yo being a shock location. Since u{x,t) = —u{—x,t), it follows that 



W{s,^+,^.,x,t) = W{s,C-,C+,x,t). 
V{s,^,t) and W{s,$,+,^-,x,t) are related by (recall that V = V+ = V- ii u{x,t) '= 



s{x, yo, t) = u{yo + x/2, t) - u{yo - x/2, t), 



and let W{s,^+,^-,x,t) be the PDF of 



{s{x,yo,t),^{yo + x/2,t),Ciyo - x/2,t)) 



-u{-x,t)) 



(3.23) 




Thus, 



(3.24) 




We have 



Theorem 3.8 



W satisfies 



{pW)t 



-psW^ + 2p{Bo - B{x))Wss 
+pi?i(I^5^5^ + W^_^_) + 2pB^{x)Wi:^^_ 

-pB2{x){Wsi+ + Wsi_) + ?1 - ?2 + J, 



(3.25) 
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with Bi{x) = —Bxx{x), B2{x) = B^ix). <ji(s,^+,^_,a;,t) is defined such that 

a;, t)dsd^-d^+dzdt, 

gives the average number of shock creation points in [z, z + dz) x [t,t + dt) with 

s{x,yi,ti) £ [s,s + ds), 

ayi-x/2,ti) G [C-,C-+d^.), 

conditional on {yi,ti) G {[z,z + dz) x [t,t + dt)) being a point of shock creation (because of 
the statistical homogeneity, z is a dummy variable). (;2is,£,+ ,(,-,x,t) is defined such that 

<i2{s, X, t)dsd(,-d^-\-dzdt, 

gives the average number of shock collision points in [z, z + dz) x [t,t + dt) with 

s{x,y2,t2) G [s,s + ds), 
i{y2 + x/2,t2) G + 
ay2-xl2M) e [i-,i- + di.), 

conditional on (1/25^2) S {[z,z + dz) x [t, t + dt)) being a point of shock collision. Fi- 
nally, J(s, ^_|_, accounts for the possibility of having another shock in between 
[yj — x/2, yj + x/2] and satisfies 

J{s,^+,C-,x,t) = 0{x). 

At statistical steady state, the definitions for iji, ^2 simplify. Indeed, in the limit as 
t — >■ +00, we have 

'^i{s,C+,(-,x,t) aiSi{s,(+,C-,x), 

where ai is the space-time number density of shock creation points, 5i(s, x) is the 
PDF of 

{s{x, yi,ti),^{yi + x/2, ti), - x/2, ti)), 

conditional on a shock being created at and 

'^2{s,C+,^-,x,t) a2S2{s,^+,^-,x), 

where a2 is the space-time number density of shock collision points, 5*2(5, ^+,^_,x) is the 
PDF of 

{s{x,y2,t2),^{y2 + x/2,t2),^{y2 - x/2,t2)), 



STATISTICAL THEORY FOR BURGERS EQUATION 
conditional on a two shocks colliding at {y^it^)- 
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Remark on the strategy for the proof of Theorem 3.8. Ideally to prove Theorem 3.8 
we should follow the strategy in Sec. 3 for the derivation of the equation for Q. Let 
X(ui,xi, . . . ,U6,X6,i) be the PDF of 

{u{yQ + xi,t), . . . ,u{yQ + XQ,t), 

conditional on yo being a shock location. Knowing the equation for X, one can easily 
derive an equation for the conditional PDF of 

{u{yo + x/2, t), u{yo + x/2, t), r]{yo + x/2, y, t),ri{yo - x/2, y, t) 

where rj(x,z,t) = {u(x + z,t) — u{z,t))/z. Letting z ^ 0, one derives an equation for 
W. Clearly, this derivation is rather tedious and, as we now show, unnecessary for our 
purpose. 

Recall that 



pX{ui,xi, . . .,ue,xe,t) = 




The average under the integral is the characteristic function associated with X, and an 
equation for this quantity can be derived using the equation for u{z + Xp,t) 

du = - [u]^Uxpdt + dW{z + Xp,t), p=l,...,6. 

Instead of reproducing these straightforward calculations, we note simply that for homo- 
geneous situations the resulting equation for X will contain terms proportional to 

P2{xp,t) = (^5{xp + z-yk)5{z-yj)Y 

These terms account for the probability of having another shock, say y\ , between yo and 
yo + Xp'. they are the origin of J in (3.25). Note also that technically, the p2{xp,t)''s arise 
because of the average [u]^ in the equation for u{z + Xp,t). Now, the key point is to note 
that 

pixp,t) = 0{xp). 



As a direct result. 



J{s,C+,C-,x,t) = 0{x). 
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We are eventually interested in the limit as x ^ of W. As argued in Sec. 3.6, in 
this limit the 0{x) terms in (3.25) are negligible. Thus, we will not dwell on obtaining an 
explicit expression for J. Instead, we will derive (3.25) using {W{x,t) = Wx{x,t)) 



du = —uux±dt + dW{z + x±,t), 

= -Hx± + ^^)dt + dW{z + x±,t), 



as if no shock are present between z and z + x±. The error we are making are accounted 
for by the term J. 

Proof of Theorem 3.8: Define 

^—^X^u(z+x+,t)—^X-u{z+X-,t)—i|J,+^{z+X-^.,t)—^|J,-^{z+X-,t) 



Then 



(27r) 

^('^' a:/2, -x/2, z, t)S{z - yj)). 



We now derive equations for (0) = (^j05{z — yj)), W. We use the following rules 
from Ito calculus 

dW{x,t)dW{y,t) = 2B{x-y)dt, 

dW{x,t)dW{y,t) = 2Bi{x - y)dt 
dW{x,t)dW{y,t) = -2B{x-y)dt, 

where Bi{x) = —B^xix), B2{x) = Bx{x). Thus, from (3.26), we obtain 

dO = i(A+ti+M+^^ + X-U-U-^ )0dt 

- {XlBo + X'iBo + 2A+A_B(.T+ - x^))6dt 

+ i{lJ^+u+i+x+ + ^^+&r + f^-'^-^-x- + fJ'~^-)Odt 

- {ix\Bi + i?_Bi + 2^+^_5i(x+ - x-))edt 

- (A_/i+ — X+ii-)B2{x+ — x-)Odt 

- iX+edW{z + x+, t) - iX-edW{z + x-,t) 

- iiJLjf.9dW{z + t) — ijJL-9dW{z + x_, t), 

where u± = u{z + x±,t), ^± = (^{z + x±,t). Similarly, using dyj/dt = u{yj,t), we get 

+ Y,6{z- yk)6{t - tk)dt - ^ <5(z - yi)6{t - ti)dt, 
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where S^{z) = dd{z)/dz, the (yfe,tfe)'s are the points of shock creation, the {yi,tiys are the 
points of shock colHsions. Using 

0S^{z - Vj) = {6S{z - yj))z - Bx+5{z - yj) - 0x_6{z - yj), 

and noting that = by statistical homogeneity, it follows that 

{@)t = iX+{u+u+^^@) + iX-{u-u^^_e) 

- {XIBq + XIBq + 2A+A_B(x+ - x_))(e) 

+ + C+)e) + + e)e) 

- {filBi + filBi + 2fi+fi^Bi{x+ - x_))(e) 

- (A_/x+ - X+fi-)B2{x+ - x-){Q) 

j 

where Si and S2 account respectively for shock creation and collision events. These are 
given by 

= e-'^+''+-'^-"^'+^+-'''-^-S{z - yk)5{t - tk) 

^2(A+,A_, //+,//_, 
I 

To average the convective terms we use: 

ix±{u±u±^^e) + ifi±{iu±^±^^ + ^|)e) 

For the term involving u{yj,t) = {u+{yj,t) + u-.{yj,t))/2 we note that 

u+{yj,t)9x+ = u{yj + x+, t)ex+ -x+ J^dp ^{yj + Px+,t)ex+ 
= {u{yj + x+, t)9)x+ - i{yj + x+, t)e 



-x+ I dp i{yj + Px+,t)9x. 
Jo 

= iOx+\+ - iOiJ,+ -x+ I dp ^{yj + Px+,t)ex+. 
Jo 
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A similar expression holds for u-{yj,t)6x_- Also 

u+{yj,t)6x_ = {u{yj + x+,t)ej)x_ -x+ [ dp £,{yj + /3x+,t)dx- 

= idx-X+ -x+ dp ^{yj + Px+,t)ex_, 
Jo 

and a similar expression holds for u-{yj,t)9x_^_- Thus 



2(y2Hyj,t)i0x++0x.)Siz-yj) 

3 

= i{Q)x+\+ + i{@)x-x+ + i{Q)x+X- + i{Q)x-X- - - - R 

where 

i?(A+, A_,/i+,/i_,x+,a;_) 
= ^+(E/^/^ i{yj+Px+,t){ex^ + ex_)5{z - yji 

+ ^-(e/^/^ ayj+Px-,t){0x+ + 0x-)S{z - yj)^ 

Combining the above expressions leads to the following equation for (G): 

(0)t = -^((0)x+A+ + (e)._A_ - (e).^A_ - (e)._Aj 

-{XIBq + XLBo + 2X+\_B{x+ - x-)){@) 

-{filBi + nlBi + 2/x+/x_Si(x+ - a;_))(e) 
-{X-H+ - \+H-)B2{x+ - ar_)(e) + Si - E2 - R. 

To obtain an equation for we note the following remarkable property of R: 

Lemma 3.9 

i?(A+,A_,/x+,/x_,x/2,-a;/2) =0. 

Proof: To see this, write 

i?(A+, X-,ii+,n-,x/2, -x/2) 
= \ 1™, ^(E ll'^i^ i^(yj + t) - ^yj - /3a^/2, t))e6{z - y,)y 
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where 6 = Oj{X+, A_,/i+, |U_,x + x/2, x — x/2, z,t). We claim that 

^ = (E/V iaVj + Px/2,t) - - Px/2,t))e6{z - y,)^ = 0. 

Indeed the symmetry u{x, t) '= —u{—x, t) requires that A be invariant under the trans- 
formation 

z -z, X ^ X, X -X, yj -yj, A± -A=p, //± fi^. 

On the other hand one checks explicitly that A — > —A under the same transformation. 
Hence A = 0, R = 0. I 

We now continue with the proof of Lemma 3.8. Combining the above expressions, on 
the subset A+ = A, A_ = —A, x+ = x/2, x- = —x/2, (6) satisfies 

= -i(G),A - 2X\Bo - B{x)){e) 

-{filBi + nlBi + 2/x+/x_Si(x))(G) 
+{XH+ + Xn-.)B2{x){Q) + E2, 

where the S\, S2 are evaluated at A+ = A, A_ = —A, x+ = x/2,x^ = —x/2. Going to the 
variables (s,^+,^_) we obtain (3.25). 

I 

3.6 The exponent 7/2 

In this section we will derive the following result: 

For large negative ^, F behaves as 

(3.27) F{^) ~ C|er^/^ as^^ -00. 



A direct consequence of (3.27) is 
(3.28) Q^iO - 



C_|Cr^/2 as €^-00, 

C+^e-^ as C ^ +00, 



The argument for (3.27) is based on the following result 
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Consider 

x-^Si{sx^/^, ^-x~'^/^,x), 
the PDF of the resettled variables 

{s{yi,x, ^l)x-l/^ + x/2, ti)x^/^^-{yi - x/2, ti)x^/^), 

conditional on a shock being created at {yi,ti). In the limit as x ^ 0, 

where P(-) is a PDF supported on (— oo,0]. 

(3.29) shows that, in the original variables, >Si is asymptotically 

^' ' ~x-^/3p(sxV3)<5(^+-sx-V3M(C+-^-). 
We first derive (3.27), then (3.29). 
Derivation of (3.27). Recall that 



F{i)=p lim / dsd^' Woo{s,i',i,x), 



where is the statistical steady state value of W. Thus, evaluating F amounts to 
evaluating W which we will do by analyzing (3.25). We note first that this equation 
describe the process of shock creation, motion, then collision. Since collision only occurs 
if shocks are present, it is natural to represent the effect of collision as proportional to the 
density of shocks in the systems, i.e. we write 



?2 = pg{s,^+,^-,x)W, 



for some function g{s, ^_|_, x) which is assumed to be smooth in s, x. This amounts 
to assuming that the characteristics of the shocks around the collision points are not very 
different from the characteristics of the shocks away from the collision points. For instance, 
if they were identical, we would have (^2 = (T2W and, hence, g = (J2I P- Next, since we are 
interested in W evaluated at x = 0, we neglect the terms 

B2{x) = 0{x), 2{Bo - B{x)) = 0{x^) J = 0{x), 

in (3.25). Finally, since wc arc interested in the limit as ^± — > —00, we neglect the forcing 
terms in proportional to Bi or Bi(x). 
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Under these approximations, (3.25) reduces to 

{pW)t = -psw, + + ^-)w + pie+w)^^ + piew)^_ - pgw + ^1. 

Assuming no shocks are present at the initial time, the equation must be solved with the 
initial condition pW{s,^+,^-,x,0) = 0. The solution is 



pW{s,^+,C-,x,t) = fdr ((l-e+T)(l-e-T))-^/2 



X exp (- /V , 3^^, - - 

The statistical steady state solution is obtained in the limit as i — > oo of this expression. 
Using 

lim <Ji(s,^+,^_,x,i) = aiS'i(s,^+,C_,a;), 

t— >oo 

we obtain 



roo 

pW^{s,^+,i^,x) = cTi / dr ((1 -e+r)(l -C_r))-^/2 

Jo 



exp (- f^dr' g ^^^^ - " -'^ 



At X = 0, using (3.30) for S\, we get 

pw^oo(s,e+,c-,o) 

= ai/ dr ((l-C+r)(l-C_r))-V2 

X exp(-^V ,,(,,. ^-i±-;,^^J^,-r'.)) 



Since 



1/3/ Vl-e+r 3t; Vl-C+r 



Mr^-r^) = (i-c.r)^^(^.-e-), 
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we have Woo{s,^+,^-,0) oc 6{^+ — ^_). Using the relation (3.23) between W and V, this 
implies that Woo(s, C+) 0) = ~^-)) ^^'^ leads to 

pV^{s,0 = a, fdr (1 -er)-3(|.|r)-V3p(_f!(!)<5 + ^) 

X exp (- /V , -r'.)) . 

To perform the integration over r we use 

Since we are considering ^ ^ — 1, the exponential factor evaluated at r = —1/2^ is 

exp (- f'^' , = 1 + O(r'). 

This means that shock collision events make no contribution to leading order, leaving us 
with 

py^(s,0 = Cai|.|-^/3|^|-5/3p(_(2,2|^|)i/3^^ 
where C = 24/3/3. Hence, 

(3.31) F{0 = -Ca, r ds .2/3|^|-5/3p (-(2.2|^|)i/3\ = -C\C\-'/^, 

where C = 2-VVi dh \h\^l'^P{h). 



DerivaUon of (3.29). We use local analysis around the shock creation points [20]. 
Consider a shock created at yi at time ti with velocity ui = u{yi,ti). Assuming x is an 
analytical function of u, we have locally 

(3.32) X = a{u{yi + x/2, h) - mf + 0{u{yi + x/2, h) - m)^), 

where a < is a random quantity. Setting a = (2/6)^ gives 

uiy, + x/2, h) = m + ^xV3 + 0(;^2/3)_ 

Hence 

s{yi,x, ti) = u{yi + x/2, ti) - u{yi - x/2, ti) 
= bx'/' + 0{x'/% 
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and 

(3.34) ^{y, + x/2, ti) = + 0(x-i/3). 

Note that these formulae are only valid if there is no other shock in [yi — x/2, yi + x/2]. 
Since the probability of having another shock in [yi — x/2, yi + x/2] is at most 0{x), the 
errors we incur by using (3.33)- (3.34) are of higher order. 
Recall that 

5i(^,e+,C-,x) = / '^'^(2%3^" e^^-+^^+^++^/--^-^^(A, x), 

where 

(Tiri(A,/x+,/i_,a;) 

^g-«As(2/fc,a;,tfc)-i//+C(j/fc+a;/2,tfc)-i//_^(2/fc-a;/2,tfc)^^^ _ y^^^^f _ 
k ' 

O is the characteristic function associated with Si. Similarly 

J^(Aa;-V3,/x+x2/3,/x_x2/3,x), 
is the characteristic function associated with the rescaled PDF 

xSi{sx^'^ ,i+x-^l^ ,i-x-^'^ ,x). 

We evaluate Q.{\x~^/^ ^Ji+x"^/^ ^Ji-x"^/^ ,x) in the limit as a; ^ using (3.33) and (3.34) for 
s{yi,x,ti), ^{yi + x/2,ti). This gives 

k ' 

In the limit as a; — > 0, b, (yk,tk) are the only random quantities to be averaged over. 
Furthermore b is statistically independent of {yk,tk) because of statistical homogeneity 
and stationarity. Let P{b) be the PDF of b. Then 

(Ti lim 0(Ax~^/^, /i+x^/^, x) 



^6{z-yk)S{t-tk)) f db p(^,)e-A^,-^(A++/i_W3 
k ' 

= ai [ db p(^b)e-'~^^-'^P++f- 

J —oo 



_)6/3 



Direct evaluation of this expression gives (3.29) in the variables (s,^_|_,^_). 
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3.7 Connection with the geometric picture 

Here we compute directly the contribution to F in the neighborhood of shock creation. 
This is a reformulation of the argument presented in [13] in terms of quantities defined in 
the present paper. Assume a shock is created at time t = 0, position x = yi, and with 
velocity u = u\. Then locally (compare (3.32)) 

X = ui + {u — uijt + a{u — uif + 0{{u — ui)^t), 

where a < is a random quantity. For the purpose of comparison with (3.31), it is useful 
to set a = (2/6)^. Since for t ^ 1 to leading order the shock is located at a; = yi, to 
leading order u-{yi,t), tt+(yi,i) are solution of = (u — ui)t — {2{u — ui)/b)^. Thus 

<'±(!/i,*)=«i=f(J^)'" + 0(*), 
(3.35) .fe,t) = -(l|!)-'%o(*) 

Similarly, to leading order ^_(yi,t), ^^{yi,t) are solution of 1 = — 3(2/6)^(u± — ui)'^^. 
Thus 

(3-36) ^^(y,,i) = _l + o(i). 



Recall that from (3.24) (using ^_|_(x,t) '= S,-{—x,t)) 



FiO = J ^e^^^{j:s{z,t)e-^^^^(^''^S{z-yj, 
'j24z,m-^+{z,mz-yj)), 



Under the assumption of ergodicity with respect to time-translation, F{^) can be evaluated 
from 



L,T-^+oo 2L1 Jo J-L J 
i,T-»+oo 2L1 Jo ~^ 



where N(t) is number of shocks in [—L, L] at time t. The contribution to F near shock 
creation points, say Fi, can be evaluated for large negative ^ using (3.35), (3.36) for 
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s{yi,t), This gives in the hmit as ^ — > — cx) 

where C = 2-^/'^ai J^^db\b\^/'^P{b). Comparing with (3.31), we conclude that Fi{^) = 
F{^) to leading order. 

4 Conclusions 

To recapitulate the highlights of this paper, by writing down and working with the master 
equations in the inviscid limit, we have shown that the scaling of the structure functions is 
related to the shocks which are the singular structure in the limiting flow. The scaling of 
the PDFs, on the other hand, is related to the shock creation and collision points, which 
are singularities on the singular structures. 

The present paper provides a framework within which various statistical quantities 
of the stochastic Burgers equation can be calculated using self-consistent asymptotics 
without making closure assumptions. The main examples used here are the asymptotic 
behavior of structure functions and the PDF of the velocity gradient. It seems likely that 
other statistical quantities, such as the tails of the velocity PDF and the PDF for velocity 
difference, can also be analyzed in the present framework by exploiting further the source 
terms in (2.6), (2.23). 

Appendix: Master equations for the viscous case 

In this Appendix we list results for the PDFs in the viscous case. The master equations 
below were previously derived, e.g., in [23, 24, 33]. 

Let {u,^,x,t) be the PDF of {u{x,t),^{x,t)) for solutions of (1.1). First we have 

Lemma A.l satisfies 

-i^{{u,,\u,opnu-''ii^xx\u,opni, 

where Bq = -B(O), Bi = —Bxx{^), {^u^^) denotes the conditional average on u and ^. 

(A.l) is unclosed since the form of {uxx\y',i) and {ixx\u,5,) entering the viscous terms 
is unknown. Most work has resorted to various closure assumptions. Our main goal has 
been to find ways to extract information from the master equations such as (A.l), without 
making any closure assumption. Note that from the identity 

Pxx = -{{^xx\n,op'')u - {{^xx\u,opn^ 

+{{ui\u,opiuu + 2{{ux^x\u,opnu^ + {{ex\y',opi^^, 
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the viscous term in (A.l) can also be written as (using {u^\u,0 = {ux^x\u,0 — 

Thus, viscous effects give rise to anti- diffusion terms since > and (^^|it, C) ^ 0- this 
is natural since viscosity tends to shrink the distribution P*^ towards the origin. 

Since ^ = Ux, we have {a{u))x = (a„('u)0 for all smooth and compactly supported 
functions a(-). This is expressed as 

Lemma A. 2 The consistency relation 

(A.2) Jd^P^ + /daP: = 0, 

holds for all time for the solution of (A.l) if it holds initially. 

Lemma A.2 can be proven upon noting that A = f d^P^ + f d^ satisfies At = 

-uAa: + BqA 

uui equation that can obtained by integration of (A.l). Since ^4 = 
initially, it is zero for all time. In the statistically homogeneous case, (A.2) reduces to 
/ d^^P^ = (or equivalently {a{u))x = {0'u{u)$.) = 0). (A.2) also ensures that this 
equation preserves the normalization of P*^. In fact we have 

Corollary A. 3 The solution of (A.l) satisfies 

(A.3) 4 / dud^ P'' = 0. 

dt J 

i.e. fdudCP'' = Jdud(Po. 

(A.3) follows immediately from integrating (A.l): 

^ J dud$, P^ = —{ux) + {Cj + boundary terms = 0. 

The first two terms at the right hand-side cancel because of (A.3); the boundary terms 
vanish because for finite v, P" decays faster than algebraically in ^ as |^| — +oo. 
Consider the reduced distributions 



Q''i^,x,t)= J duP''{u,^,x,t), 



R''{u,x,t) = / d^ P''{u,^,x,t). 
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(A.2), written after integration over u as 

r du' Rl{u',x,t)+ f d^^P'' = 0, 

J — CO J 

can be used to derive from (A.l) an equation for R^: 

4^ 1^ = -uK - /"^ du' Rl{u\ X, t) 

+BoK^-u{{u,.\u)Rn^. 

For statistically homogeneous situation, = in (A.l). Then, using / d^^^P'^ = and 
(A.l) leads to the following equations for and Q^: 

(A.5) = So<„ - lyiiu^MRlu- 



(A.6) = ^Q-" + ieQ^k + BiQii - mxxlOQlc 
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